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Abstract 

We consider the problem of the nature and possible types of spacetime singular- 
ities that can form during the evolution of FRW universes in general relativity. We 
show that by using, in addition to the Hubble expansion rate and the scale factor, 
the Bel-Robinson energy of these universes we can consistently distinguish between 
the possible different types of singularities and arrive at a complete classification of 
the singularities that can occur in the isotropic case. We also use the Bel-Robinson 
energy to prove that known behaviours of exact flat isotropic universes with given 
singularities are generic in the sense that they hold true in every type of spatial 
geometry. 
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1 Introduction 



There has recently been a resurgence of interest in tlie old question of the existence 
and nature of spacetime singularities, especially in a cosmological context. Tradition- 
ally one is interested in formulating criteria, of a geometric nature, to test under what 
circumstances such singularities, in the form of geodesic incompleteness, will be formed 
during the evolution in general relativity and in other metric theories of gravity, cf. p. 
These are usually translated into sufficient conditions to be satisfied by the matter fields 
present and are very plausible. On the other hand, one can formulate equally plausi- 
ble and generic geometric criteria for the long time existence of geodesically complete, 
generic spacetimes in general relativity and also other theories of gravity, cf. 0. Such 
criteria assume a globally hyperbolic spacetime in the so-called sliced form (cf. 
require the space gradient of the lapse function as well as the extrinsic curvature not to 
grow without bound, and are also very plausible. 

The lesson to be drawn from this state of affairs is that in general relativity (and also 
in other metric gravity theories), singular and complete spacetimes are equally generic 
in a sense. The real question is what do we mean when we say that a relativistic 
model develops a spacetime singularity during its evolution. In other words, what are 
the possible spacetime singularities which are allowed in gravity theories? These may 
be singularities in the form of geodesic incompleteness developing in the course of the 
evolution but also others which are more subtle and which are also dynamical and present 
themselves, for instance, in some higher derivatives of the metric functions spoiling the 
smoothness of global solutions to the field equations. 

Obviously a recognition and complete analysis of such a program cannot be accom- 
plished in the short run and requires complete examination of a number of different 
factors controlling the resulting behaviour. For instance, one needs to have control of 
each possible behaviour of the different families of relativistic geometry coupled to matter 
fields in general relativity and other theories of gravity to charter the possible singular- 
ity formation. In this sense, unraveling the nature and kinds of possible singularities in 
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the simplest kinds of geometry becomes equally important as examining this problem 
in the most general solutions to the Einstein equations. In fact, an examination of the 
literature reveals that the more general the spacetime geometry considered (and thus 
the more complex the system of equations to be examined) the simplest is the type of 
singularities allowed to be examined. By starting with a simple cosmological spacetime 
we allow all possible types of singularities to come to surface and be analyzed. 

In jlj we derived necessary conditions for the existence of finite time singularities 
in globally and regularly hyperbolic isotropic universes, and provided first evidence for 
their nature based entirely on the behaviour of the Hubble parameter H = a/ a. This 
result may be summarized as follows: 

Theorem 1.1 Necessary conditions for the existence of finite time singularities in glob- 
ally hyperbolic, regularly hyperbolic FRW universes are: 

51 For each finite t, H is non-integrable on [ti,t], or 

52 H blows up in a finite future time, or 

53 H is defined and integrable for only a finite proper time interval. 

Condition Si may describe different types of singularities. For instance, it describes 
a big bang type of singularity, when H blows up at ti since then it is not integrable on 
any interval of the form [ti,t], t > ti (regular hyperbolicity is violated in this case, but 
the scale factor is bounded from above). However, under 5*1 we can have other types of 
singularities: Since H{t) is integrable on an interval [ti,t], if H{t) is defined on [ti,t], 
continuous on (ti, t) and the limits lim^_^j+ H{t) and lim^_^f- H{t) exist, the violation of 
any of these conditions leads to a singularity that is not of the big bang type discussed 
previously. 

Condition 5*2 describes a future blow up singularity and condition S3 may lead to 
a sudden singularity (where H remains finite), but for this to be a genuine type of 
singularity, in the sense of geodesic incompleteness, one needs to demonstrate that the 
metric is non-extendible to a larger interval. 
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Note that these three conditions are not overlapping, for example Si is not implied 
by 5*2 for if H blows up at some finite time tg after ti, then it may still be integrable on 

[tut], ti<t < ts. 

There are many examples in the literature of singularities belonging to the types 
predicted by the above theorem; they were analyzed in ^ in detail. They usually 
describe flat isotropic universes with various components of matter. For example, those 
models containing phantom dark energy 5J, or a perfect fluid and a scalar field [Hj, 
exhibit future finite time blow up singularities and therefore fall in the 5*2 category. The 
standard big bang singularities as well as the sudden singularity of [7j fall in the Si 
category. Other sudden singularities (which do not have a blow up of at t = 0, see for 
example [Hj) and the inflation model of |2j both fall in the S^ category. 

2 Classification 

As discussed in PJ, Theorem fll.ip describes possible time singularities that are met in 
FRW universes having a Hubble parameter that behaves like 5*1 or 5*2 or 6*3. Although 
such a classification is a first step to clearly distinguish between the various types of 
singularities that can occur in such universes, it does not bring out some of the essential 
features of the dynamics that differ from singularity to singularity. For instance, con- 
dition 5*2 includes both a collapse singularity, where a ^ as t — > tg, and a blow up 
singularity where a — > 00 as t ^ t^. Such a degeneration is unwanted in any classification 
of the singularities that can occur in the model universes in question. 

It is therefore necessary to refine this classification by considering also the behavior 
of the scale factor. Another aspect of the problem that must be taken into account 
is the relative behaviour of the various matter components as we approach the time 
singularities. 

In what follows we present a summary of our most recent work on the classification of 
possible singularities which occur in isotropic model universes. Details will be published 
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elsewhere, cf. ^U]. Our classification is based on the introduction of an invariant geomet- 
ric quantity, the Bel-Robinson energy, which takes into account precisely those features 
of the problem, related to the matter contribution, in which models still differ near the 
time singularity while having similar behaviours of a and H . In this way, we arrive at a 
complete classification of the possible cosmological singularities in the isotropic case. 

The Bel-Robinson energy is a kind of energy of the gravitational field projected in a 
sense to a slice in spacetime. It is used in jTT], [12] to prove global existence results in the 
case of an asymptotically fiat and cosmological spacetimes respectively, and is defined 
as follows. Consider a sliced spacetime with metric 

in+l)g^_f^2^Qy ^ ^0 = rft, e' = dx' + P'dt, (2.1) 

where N = N{t,x^) is the lapse function and l3^{t,x^) is the shift function, and the 
2-covariant spatial electric and magnetic tensors 

pi _ pO 
— rtjQj , 

n — —n n jDhklm 



hk 
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where rjijk is the volume element of the space metric g. The Bel-Robinson energy is given 
by 

Bit) = \! [\E\' + \D\' + \B\' + \H\') (2.2) 

JMt 

where by |Xp = g^^ g^^XikXji we denote the spatial norm of the 2-covariant tensor X. 
In the following, we exclusively use an FRW universe filled with various forms of matter 
with metric given by 

ds^ = -dt^ + a^{t)d(T'^, (2.3) 

where da"^ denotes the usual time-independent metric on the 3-slices of constant curvature 
k. For this spacetime, we find that the norms of the magnetic parts, \H\, \B\, are 
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identically zero while \E\ and \D\, the norms of the electric parts, reduce to the forms 

|E|2 = 3 (a/a)' and \D\^ = 3 {{d/af + k/a'Y- (2.4) 

Therefore the Bel-Robinson energy becomes 

B{t) = ^ {\E\' + \D\') , (2.5) 

where C is the constant volume of (or in in the case of a non-closed space) the 3- 
dimensional slice at time t. 

We can now classify the possible types of singularities that are formed in an FRW 
geometry during its cosmic evolution and enumerate the possible types that result from 
the different combinations of the three main functions in the problem, namely, the scale 
factor a, the Hubble expansion rate H and the Bel Robinson energy B. If we suppose 
that the model has a finite time singularity at t = tg, then the possible behaviours of the 
functions in the triplet [H, a, (|ii^|, |D|)) in accordance with Theorem are as follows: 

H non-integrable on [ti,t] for every t > ti 

5*2 H ^ oo at ts > ti 

5*3 H otherwise pathological 

iVi a^O 

N2 a ^ as ^ 

N3 a ^ 00 

Bi \E\ — s> 00, \D\ 00 
B2 \E\ < 00, \D\ 00 
i?3 \E\ 00, \D\ < 00 
B4 \E\ < 00, \D\ < 00. 
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There are a few types that cannot occur. For instance, we cannot have an (S'2, A^2, -B3) 
singularity because that would imply having a < 00 (A^2) and H — > 00 {S2), while 
3 {{a/aY + kjd^Y < 00 (-B3), at tg which is impossible since 00 at ts {k arbitrary). 

A complete list of impossible singularities is given by {Si, Nj, Bi), where the indices in 
the case of a A; = 0, ±1 universe take the values i = 1,2, j = 1,2,3, / = 3,4, whereas 
in the case of a A; = —1 universe the indices take the values i = 1,2, j = 2,3, I = 3,4 
(here by Si we denote for simplicity only the big bang case in the Si category). We thus 
see that some singularities which are impossible for a flat or a closed universe become 
possible for a hyperbolic universe. Consider for example the triplet {S2, Ni, B^) which 
means having — > 00, a — > and 

\D\^ ^3 {{d/ af + k/a^)^ < 00 

at tg. This behaviour is vahd only for some cases of a hyperbohc universe. 

All other types of finite time singularities can in principle be formed during the 
evolution of FRW, matter-filled models, in general relativity or other metric theories of 
gravity. 

It is interesting to note that all the standard dust or radiation-filled big-bang sin- 
gularities fall under the strongest singularity type, namely, the type {Si, Ni, Bi). For 
example, in a flat universe flUed with dust, at t — we have 

a{t) oc t^/^^O, (iVi), (2.6) 
H oc r^^oo, {Si), (2.7) 
\E\^ = 3/AH^^oo, = 3i/^ ^ 00, {Bi). (2.8) 

Note that the classification is organized in such a way that the character of the singu- 
larities becomes milder as the indices of S, N and B increase. Milder singularities are 
thus met as one proceeds down the list. 

In fluid-filled models, the various behaviours of the Bel-Robinson energy density can 
be related to four conditions imposed on the density and pressure of the cosmological 
fluid: 
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Bi fi ^ CO and |/i + 3p\ ^ oo 
B-2 -v^ /i — >• oo and |/i + 3p| < cxo 

i?3 -vv- /i < oo and |/i + 3p| — cxd < oo and |p| ^ oo 
^4 -x^ < oc and |/i + 3p| < oo <^=^ /i < oo and |p| < oo. 

Of course we can translate these conditions to asymptotic behaviours in terms of a, H, 
depending on the value of k, for example, 

1. li k = 0, fi < oo => < oo, a arbitrary 

2. If A; = 1, /X < oo ^ iJ^ < oo and a 7^ 

3. If = -1, /i < 00 ^ if^ - < 00. 

As an example, we consider the sudden singularity introduced in ||7j. This has a finite 
a (condition A^2), finite H (condition 5*3), finite ^ but a divergent p (condition B^) at 
tg. As another example, consider the flat FRW model containing dust and a scalar field 
studied in [Hj. The scale factor collapses at both an initial (big bang) and a final time (big 
crunch). The Hubble parameter and a/ a both blow up at the times of the big bang and 
big crunch (cf. jlj) leading to an {Si, Ni, Bi) big bang singularity and an Ni, Bi) 
big crunch singularity, respectively. 

3 Generic results and examples 

In this Section we provide necessary and sufficient conditions for the occurrence of some 
of the singularities introduced above. These conditions are motivated from studies of 
cosmological models described by exact solutions in the recent literature. By exact 
solutions we mean those in which all arbitrary constants have been given fixed values. 
We expect the proofs of these results to be all quite straightforward, for we have now 
already identified the type of singularity that we are looking for in accordance with our 



7 



classification. Proving such results without this knowledge would have been a problem 
of quite a different order. 

The usefulness of the results proved below lies in that they answer the question of 
whether the behaviours met in known models described by exact solutions (which as a 
rule have a flat spatial metric {k = 0)) continue to valid in universes having all possible 
values of k described by solutions which are more general than exact in the sense that 
some or all of the arbitrary constants still remain arbitrary. See P| for more results of 
this type. 

Theorem 3.1 Necessary and sufficient conditions for an {S2, N^, Bi) singularity occur- 
ring at the finite future time tg in an isotropic universe filled with a fluid with equation 
of state p = wfi, are that w < ~1 and \p\ -^00 at tg- 

Proof. 

Substituting the equation of state p = wfi in the continuity equation fi + 3H{fi + p) = 
0, we have 

/X oc a~^^'"+^\ (3.1) 
and so if w < — 1 and p blows up at tg, a also blows up at tg. Since 



we conclude that at tg, H, a, \D\ and \E\ are divergent. 

Conversely, assuming an {S2, N^, Bi) singularity at tg in an FRW universe with the 
equation of state p = wfi, we have from the (Bi) hypothesis that fi —* 00 at tg and so p 
also blows up at tg. Since a is divergent as well, we see from 1)3.11) that w < —1. ■ 

As an example, consider an exact solution which describes a fiat, isotropic phantom 
dark energy filled universe, studied in ^ given by 

2 

a = [a,'^'^-y^ + Mi±^(t _ toi"^' , (3.3) 
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where A is a constant. From ()3.1|) we see that the scale factor, and consequently if, 
blows up at the finite time 

2 



ts — to + 



3VA{\w\ - l)al 



3(|t«|-l)/2- 



Then it follows that lE'p = + 3w)^ and |Dp = 3H^ also blow up at ts- Therefore 

in this model the finite time singularity is of type (5*2, N3, Bi). 

The following result says that the strongest big bang type singularities that can occur 
in universes with a massless scalar field are due to the kinetic term. 

Theorem 3.2 A necessary and sufficient condition for an [Si, Ni, Bi) singularity at ti 
in an isotropic universe with a massless scalar field is that ^ 00 atti. 

Proof. From the continuity equation, (l) + 3H(J) = 0, we have (f) oc a~^, and if (p ^ 00 
then a 0. Since 

^^-f-^-- (3,4) 

H becomes unbounded at ti. In addition, since 

3 ~ 12 

and 



D\^ = ^L = ^^oo, (3.5) 



\Ef = ^{fx + 3pY = ^^oo, (3.6) 

at ti, both \D\ and \E\ diverge there. 

Conversely, assuming an {Si, Ni,Bi) singularity at ti, we have (from Bi) that ^ 00 
and so 0^ ^ 00 at ti. ■ 

As an example of the above behaviour we can use the exact solution for a fiat isotropic 
universe given in J^]- The solution is given by 

H = -, = ±./?in-. 
3t' ^ V 3 c 
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Since a oc t^^^ we have that at t = 0, a (A^i), H = ^ oo (Si), \E\'^ = 
and |Dp = 3H'^ oo (-Bi). As it follows from jT^I, this exact solution represents the 
asymptotic behaviour of a scalar field model if 

lim e-^W{(j)) = 0. (3.7) 

0— >itoo 

4 Discussion 

We have reviewed some of the main results of our recent work [10.^ regarding the clas- 
sification of singularities in isotropic universes. The possible behaviours of the three 
functions, H, a and B exhaust the types of singularities that are possible in an isotropic 
universe. We have used known exact solutions as examples to illustrate the possible 
singularities that can be accommodated in the categories included in the classification. 
We have also shown that the Bel-Robinson energy can be used invariable as a tool to 
test the nature of singularities described by given exact solutions and decide whether 
such behaviour is generic and independent of the spatial geometry. 

It would be interesting to further investigate how this classification will change when 
we consider other types of model universes, matter fields different from the fluid ones we 
considered in this paper, or other metric theories of gravity. 

We note that the "crudest" singularity type {Si, Ni, Bi), seems to accommodate 
all standard big bang singularities known. However, the precise analytical nature of 
the other new singularity types presented here, although certainly possible, remains a 
mystery. It is only through a combination of analytic and geometrical techniques that 
more light will be shed in such interesting questions of principle. 
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